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We study sum rules for charge exchange operators of the Gamow- Teller and Fermi type and 
explore what can be learned from them about proton-neutron pairing in T — and T = 1 channels. 
We consider both energy weighted (EWSR) and non-energy weighted (NEWSR) sum rules. We 
use a schematic Hamiltonian to stress the distinctive roles of T — 1 and T = pair interactions 
and correlations in the sum rules. Numerical results for 44 Ti are presented for both schematic and 
realistic interactions. 



I. INTRODUCTION 



The pairing interaction in finite fermion systems, as are atomic nuclei, has been widely investigated since the late 
fifties [1,2]. Pairing interaction of alike nucleons has since then been taken into account in almost any microscopic 
calculation of various properties of medium and heavy nuclei. For a long time it was considered that proton-neutron 
pairing interaction may be ignored in medium and heavy nuclei due to the energy differences of proton and neutron 
Fermi levels. Despite this belief theoretical research in this field was performed [3-6]. Nowadays, at the era of heavy 
(SI . ion accelerators and of radioactive beams, allowing to study excitations of high spin states, as well as unstable nuclei 
close to the N = Z line, the problem of proton-neutron (irv) pairing is being reconsidered and is a subject of much 
"'nT , debate. 

For instance, in the case of 52 Fe isotope the data about high spin states with J > 10 can be consistently explained 
if one assumes that from there on a new band with a broken proton- neutron pair starts [7] . In particular, ttv pairing 
T = is considered to be an important source of high spin [8] by some groups. However, other groups [9] claim that 
7Ti/ pairing T — does not exist at all. 

Furthermore, T — pairing has as added interest that it is a new form of pairing, only possible in hadronic systems. 
' It is clear that, for various reasons, the subject deserves attention and any contribution which might approach the 
^ , answer concerning the existence of irv pairing in T = 1 and T = channels is welcome. 

Here we study the echo of the nv pairing of isospin T = and T = 1, seen in two sum rules for /3-decay operators. 
We consider both energy weighted (EWSR) and non-energy weighted (NEWSR) sum rules for these operators. 

The sum rule approach is a useful classical tool to get global characterizations of the spectrum and response of a 
given nucleus to particular operators [10]. The difference between (3~ and /3 + non-energy weighted sum rules {N — Z 
sum rule) was considered by several authors in connection with single (3 decay and (p, n) and (n,p) reactions [11-17]. 
In ref. [14] the charge exchange modes were studied on equal footing with the Ml mode, by comparing the results 
obtained with RPA and TDA approaches. Several approximations which might account for the missing strength of 
the Gamow- Teller (GT) states were discussed by Auerbach et al, in ref. [13]. Also effects of deformation [18,19] on 
GT strengths have been studied within a selfconsistent microscopic framework. 

In 1987, Cha considered the influence of the particle-particle channel of the two-body interaction on the (3 + strength 
[17]. Indeed, he noticed a large sensitivity of the transition amplitude on the strength of the particle-particle inter- 
action. This idea was extended by several groups to the double beta decay where all formalisms overestimated the 
transition amplitudes [20]. To correct for this drawback new approaches going beyond QRPA have been proposed. 
While the QRPA approach obeys the N — Z sum rule, any of the higher QRPA formalisms violates it. Thus, the first 
order boson expansion formalism [21] produces a sum rule which is by 20-25% larger than N — Z . The renormalized 
QRPA [22], by contrary, yields a sum rule which is smaller than TV — Z by an amount of 20%. In ref. [23] one of 
us (A.A.R.) formulated a first order boson expansion in terms of renormalized bosons and succeeded to bring the 
agreement of the sum rule with the N — Z value within 10%. 

Here we would like to see in what direction are modified the N — Z sum rules when the irv pairing in a general 
sense is taken into account. More specifically, we address the question whether one can draw useful conclusions about 
the existence of irv pairing from the EWSR and NEWSR of (3~ and (3 + strengths. 
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We accomplish this project according to the following plan. In Section 2, we treat the non-energy weighted sum 
rules and revisit Ikeda sum rule, while in Section 3 we focus on the energy weighted sum rule. In both cases the 
transition operators correspond to those of (3~ and f3 + transitions of Gamow- Teller and Fermi type. A numerical 
application for 44 Ti is presented in Section 4. Final conclusions are summarized in Section 5. 



II. NON-ENERGY WEIGHTED SUMMED STRENGTHS AND THE IKEDA SUM RULE 



The total strength for the GT operator, acting in (3 transitions and charge exchange (p, n) reactions is defined as 
follows: 



Y+ 



(2.1) 



where \I) denotes the ground state of the beta decaying nucleus, usually named as mother system, while is a 

complete set of states describing the final nucleus, which can be reached with the GT dipole operator: 



(2.2) 



Due to the completeness property of the set the summation involved in Eq. (2.1) can be performed and one 

arrives at the result: 



GT 



I). 



(2.3) 



For the sake of completeness we specify the conventions we adopt for the isospin raising operator, t + , and third 
component of the Pauli isospin matrix, t z : 



t+ \v) = |tt) ; t+ |tt) = 0; t z |tt) = |tt) ; r z \v) 



\v). 



(2.4) 



The states \n) and \v) are single particle states specified by several quantum numbers describing the motion in 
coordinate space, as well as spin and isospin degrees of freedom. The isospin quantum numbers are labeled by 7r if the 
occupying nucleon is a proton and by v if that is a neutron. Since the isospin operators do not change other quantum 
numbers but the isospin ones, we only specify the latter here. With these conventions it is clear that Bq T describes 
the total strength for the /3~ decay mode. 

Similarly, for the GT operator responsible for /3 + -decay and (n,p) reactions 



(2.5) 



the total strength is: 



B, 



GT 



Y+ ■ Y" 



I). 



(2.6) 



Although, separately, the two total rates are not easy to be calculated, their difference can be exactly evaluated with 
a minimal effort: 



B+ T - B GT = (7| ]T (at ■ Si) (trtf - tftj) \I) . 



Using commutation relations for isospin operators 



t~t + - t + t~ — -()■ T Z 



one obtains: 



B+ T -B GT = - (7| J2 (^) 2 Tf \I) - -6 (/| T z \I) =3(N-Z), 



(2.7) 



(2.8) 



(2.9) 



where TV and Z are the number of neutrons and protons, respectively, involved in the system under consideration. 
The equation obtained 
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B+ T -B GT = 3(N-Z), (2.10) 

bears the name of its founder, we refer to it as Ikeda sum rule [11]. A similar sum rule holds also for Fermi transitions 
for which the transition operators are scalars against rotations in configuration space: 

(2-11) 

i 

Repeating the procedure used above for the GT transitions, one arrives at the expression: 

B+-B~ = (N-Z). (2.12) 

We should remark the fact that there is no approximation involved in deriving these Ikeda sum rules. This 
observation infers that these sum rules are very useful in testing the accuracy of any approximate scheme for the 
description of the mother nucleus ground state as well as of the ground and excited states in the daughter nucleus. 
They are also useful experimentally to define model independent quenching of GT strengths and serve to define the 
effective qa factor in nuclei. 

In what follows we focus on how sum rules and energy weighted sum rules depend on pairing interactions. First of 
all, we notice that Y~ ■ Y+ contains a one-body and a two-body operator: 

Y- ■?+ = £ trtt + E * ' *i • (2- 13 ) 

i i=jtj 

The one-body term is simply three times the neutron number operator. This can be shown in coordinate representation 

E k4 = £ 3 (l - *0 = I ( A - 2T ^ = 3N - ( 2 - 14 ) 

i=l i=\ ^ ' 

and can also be shown in second quantized notation as follows 

A 

E *r*?" = E (ap\t-t + \a'p')at p a alpl =^o+,flav = N, (2.15) 

i—l ap.a' p' au 

where p stands for neutron (v) or proton (tt) and a for the remaining quantum numbers of the single-particle basis. 
Then, we can write 

Kt - w E * ■ & i *r*/ i J ) = 3N + ( J i E * • s i *r*/ i J ) • ( 2 - 16 ) 

Note that the sum over i,j runs over all nucleons, and the two-body term can be of the same order as the one-body 
term and of opposite sign. 
Likewise 

Y+ ■ ?~ = E 3t+tr + E ^ ■ 8 S t+tj =3Z + 5>- % t^t - (2-17) 

and 

B GT = 3Z+ (I\ E ^ ■ Sj trtf | J) . (2.18) 

Therefore, and B GT contain the same two-body term, which does not contribute to Ikeda sum rule, i.e., Ikeda 
sum rule cancels out all the two-body correlation effects in the ground state wave function. This is the reason why 
B GT depend on irv pairing but the Ikeda sum rule does not. 

The two-body term takes care of the correlations in the nuclear ground state, including Pauli correlations. Its 
evaluation requires explicit knowledge of the ground state. In particular we note that for a filled proton-neutron 
I— shell the net contribution to B GT of the one-body plus two-body terms must be zero. Thus, in shell model 
calculations Eqs. (2.16) and (2.18) can be replaced by 
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GT 



E 



(2.19) 



GT 
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(2.20) 



where iV„ and Z v denote valence nucleons outside closed I— shells and the sums over i,j run also over the valence 
nucleons only. However, in calculations based on self consistent deformed mean field like HFB and QRPA one has to 
consider the sums over all nucleons as in Eqs. (2.16) and (2.18), as protons and neutrons occupy slightly different 
orbitals and moreover with an underunity occupation probability. We note that by definition N v — Z v = N — Z. 
Similarly, for the Fermi total rate we obtain: 



and 



B + F = (I\ E *. r tf \I)=N + (I\ ]T trtf | J) , 



Bp = Z + (I\J^tZt+ \I) , 



(2.21) 



(2.22) 



where again the two-body term takes care of Pauli correlations. For the Fermi case we have that the sum of (3 + and 
(3~ strengths can also be easily calculated as 



B- 



(I\ T~T+ + T+T- \I) = 2 (I\ f 2 - T 2 |7) 



(2.23) 



Thus, if we assume that the ground state has good isospin T = \T Z \ = \N — Z\/2, we end up with the sum rule 

B+ + B F = \N - Z\ , 
which together with the Ikeda sum rule, gives the well known results 

B F =0, 



(2.24) 



B+ = N-Z 



(2.25) 



assuming N > Z, in the SU(2) limit. 

Thus, in the exact SU(2) limit, i.e., in the limit where isospin is an exact symmetry, the two-body term contributes 
subtracting a quantity equal to the maximum number of proton-neutron pairs in the A-body system. This is equal 
to the Pauli correlation. Indeed the same result is obtained in the limit of Slater determinantal wave function with 
identical proton neutron orbitals. For each tt and v occupying the same single particle state their total isospin has to 
be zero due to antisymmetrization. Because of Pauli principle the antisymmetrized product wave function contains 
the maximum possible number (Z, for N > Z) of -kv pairs coupled to T = 0, and each pair contributes with —1 to 
the two-body term. Indeed it is easy to show that for a ttv pair, the following identity holds: 



Denoting by T^ u the ttv pair isospin and taking care of the fact that the nucleon isospin is | one obtains: 



= \ (tI - 1) + \ + i (C x t v ) z . 

Now let us apply the operator t~t+, given by Eq. (2.27), on a 7w pair coupled to isospin T, (tt <g) v) TT =0 
is: 



TO/ 



\t{t + i)-\ 



2T-1 



(2.26) 

(2.27) 
The result 

(2.28) 



valid for T = 0, 1. 

Here we used the fact that the last term in Eq. (2.27) is antisymmetric and has zero expectation value in a state 
of given symmetry, 
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((77 ® IS) TQ |(C xC)J(7r®^) T0 ) =0. (2.29) 

An immediate consequence of this property is that the operator i (t w x t„) does not contribute to the matrix 
element (7| ^2 nv t~ i+ \I) if the ground state has no isospin mixing. Hence, the operator X^yj acting on the 

ground state will give a factor 1 for any nv pair coupled to T = 1 and a factor —1 for any nv pair coupled to T = 0. 
Thus, for Fermi like ^ summed strength one obtains: 

B+ = N - Z + C {itv) t , Bp = C {7Tv) t , (2.30) 

where C^"^ is a small correlation function that takes into account small isospin mixing in the ground state and that 
tends to increase with breaking of T = pairs. 

This is consistent with the fact that T = 1 wv pairs are the only ones actually contributing to the states excited by 
the Fermi operator. 

The evaluation of the two-body term for the GT strength 

B+ T = 3N + 2j2(I\^-^t^\I)- (2-31) 

7TL> 



is more involved and requires exactly solvable models to be evaluated analytically. Evaluation of B GT for valence 
nucleons in an I— shell in SU(4) and SO (5) limits can be found in Ref. [25]. 

In particular for the model discussed in Ref. [25] of M nucleon pairs in an I— shell of degeneracy 2f2, analytic 
expressions can be given of the two-body matrix element. For the case of a standard isovector spin-singlet hamiltonian 
the two-body matrix element in the ground state with isospin T = (N — Z)/2 we get 

IT\X^Z nt-i+\T\ or7 {T + l){M + T + oo)+® + oo/€ 

(i\ jL °i ■ °i l i h l J ) = ~ sz ( 2 r + 3)(n + i/2) (2 - 32) 

On the other hand, assuming a spatially symmetric correlated pair only, the contribution of the spin factor to the 
matrix element can be easily calculated and finally one gets: 



((tt <g> u) T0 \ ■ g v t n t+ | (tt ® v) T q) 
2T+1 



S{S + \)-\ 



[T(T + 1) — 1] 

, ; . for T = (S = 1) and T = 1 (S = 0) . (2.33) 
Therefore, in a simple schematic model with independent T = 1 and T = pairs, the total B^ T strength will be 

B+ T = 3N- ^W^ u)t ^ + Af^T=o | . (2.34) 
It is obvious that the f3 + strength can then be expressed as: 

B GT =3Z- {zAf^^i + aA(^)t=o I . (2.35) 

Thus, for GT transitions, increasing both T = 1 and T = number of nv pairs (AT) leads to decreasing strengths. 
The important consequence of this is that using these expressions of total strengths for Gamow- Teller and Fermi 
transitions one can gain insight on the number of pairs with T = 1 and T = 0. 

We note that for fixed number of pairs and 2T = N — Z, Bq T decreases as T increases, in agreement with the 
behavior observed in Figs. 3 and 9 of Ref. [25]. 

It is interesting to mention that these equations may also provide bounds for Gamow- Teller and Fermi transition 
total strengths. Since the total number of pairs in above equations cannot exceed the value of Z (assuming Z < N), 
from Eqs. (2.30)-(2.35) one derives: 

< B GT < 3Z, 
3{N -Z)< B+ T < 3N, 
< 6™ )r = Bp, 

(N-Z)<(N-Z)+ C (7r,y)r = B+, (2.36) 
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Comparing the second inequality of the second row in the above set of equations, with Eq. (2.31) one concludes that 
the two-body term, which is washed out by Ikeda sum rule, has a negative contribution to Bq T . Moreover neglecting 
the two body term the f3~ strength reaches its upper bound determined by the one body contribution. Thus, T = 
and T = 1 -KV pairs tend to reduce the B^ T strengths. Owing to the previous discussion on Eqs. (2.19,2.20), in 
shell model calculations the upper bounds 3Z and 3N in Eqs. (2.34) to (2.36) can be replaced by 3Z V and 3N V , 
respectively. 

III. ENERGY WEIGHTED SUM RULES 

So far, no model Hamiltonian has been used, but we implicitly assumed that there are attractive proton-neutron 
interactions in T = (S = 1) and T = 1 (S — 0) channels, favoring such pairs in the ground state. For our 
further consideration it is necessary to assume a many-body Hamiltonian H, including both one-body and two-body 
interactions. The one-body term is denoted by H and the term describing the proton-neutron interactions will be 
denoted by H p . Thus, we consider the following decomposition of H: 

H = Hq + Hp. (3.1) 

Moreover, we assume that H p is isoscalar and involves pairing in the two possible isospin channels T = and T = 1, 
which may have different strengths luq and uj\. This will allow us to discriminate between the contributions brought 
by the two terms. For pragmatic reasons it is convenient to use schematic constant interactions where only the isospin 
dependence is explicit, 

H P = ~u \ (1 - ri • r 2 ) -u>^(l + Ifi • f 2 ) . (3.2) 

Through direct calculations one proves that: 

H p (tt <g> u)\o = —ui(n ® i/) io, H p (tt <g> 7r)n = —ui(n <g> 7r)n, 
H p {v<%) u)i-t = -bj\(y ® v)x-\,H v {'k ® v) ao = -gj (tt ® v) 00 . (3.3) 

In other words, the first term of H p is a projector onto isospin T = pairs and the second term is a projector onto 
isospin T = 1 pairs. If one wants to treat the T = 1 interaction between pairs of alike nucleons with distinct strength, 
one may supplement H p with an additional term: 

ff; = -^[(i+^ 2 2 ), (3.4) 

which breaks isospin invariance. 

Here we restrict our considerations to H p which does not discriminate between different channels of the nucleon 
nucleon pairing T = 1 interaction. We assume that T = and T = 1 pairing interactions are attractive and 
consequently u>q,u>i > 0. 

In what follows we shall study the energy weighted sum rules associated to the many body Hamiltonian described 
above and f3~ transitions. For Fermi transition, EWSR is determined by the transition operator F defined by Eq. 
(2.11): 



S F+ = (EWSR) F+ ee J2(Ef Ej)\(F\F + \I)\' 



= ±(I\[F-,[H,F+]]\I) (3.5) 

The structure of the many body Hamiltonian given by Eq. (3.1) induces the following split of the energy weighted 
sum rule: 

(EWSR) F+ = 5 F+ + S F+ . (3.6) 

where the index k(= 0,p) indicates that the term corresponds to the Hamiltonian Hk- 
The only terms which survive in the first commutator of Eq. (3.5) for H = H p are: 
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.■^Etf = E ^ ^Etf + ^E^ 

J \=x,y,z y j j 

E E k A (^ a >2 fc -« 1Afc r 2 fc )^ a + (1^2)}=0. 

j \=x,y,z 



(3.7) 



where the standard notation for the antisymmetric unity tensor e^, is used. Due to this relation it is manifest that: 

[H P ,F+]=Q. (3.8) 
Therefore EWSR for the /3~ transition of Fermi type, is unchanged by the ttv pairing interaction: 



*r = o, 



(3.9) 



as it is obvious because the Fermi operator commutes with any isoscalar operator. If Hq is not isospin invariant 
(a common case when the mean field single particle states for protons and neutrons are different from each other), 
adding only an attractive T = pairing interaction to Hq will increase the number of ttv pairs coupled to T = 0. 
This decreases the total strength, i.e., B F + < N (see Eq. (2.30)), while the EWSR is left unchanged. In order that 
this picture is achieved, some of the strength is to be moved to the higher energy states. The opposite is true for 
only T — 1 -kv pairing interaction, i.e., the total strength is increased and a shift of the transition strength toward 
the low energy states is expected. Applying the hermitian conjugate operation to Eq. (3.5) one also obtains that the 
proton-neutron interaction does not affect the EWSR associated to the (3 + transition, 



Sf = 0, 



(3.10) 



and similar conclusions concerning the total strength and the migration of the strength to lower or higher energy 
states, as for the case of /3~ , hold. We stress that no migration of the strength is to be expected when H is isospin 
invariant and the ground state has isospin T = T z = (N — Z)/2 both before and after adding the pairing interaction. 

Consider now the case of GT transitions. For this case also the two terms of the many body Hamiltonian (3.1) 
determine the following decomposition of the energy weighted sum rule: 



r-GT± , cGT± 
D ' p 



The term generated by the pairing interaction has the expression: 



H Y + 



(3.11) 



(3.12) 



The double commutator involved in the above equation can be straightforwardly calculated if one uses the inter- 
mediate results: 



Y- 



T ■ T Y + 
' l ' 2 ' 1 a 



= (*1-<?2)«T+-<T+). 

= -(<?i- a 2 f (t ± ■ t 2 + t[tI) 



(3.13) 



To calculate the average of the double commutator we have to know how each of the two factors, depending on spin 
and isospin respectively, acts on a given (ttu) pair. 

((tt ® u)to\ (t • r + «) |(tt ® v) T q) = 4(T - 1), for T = 0, 1, 
((7r<g>7r)n| (t, ■ f 2 +rfr|) |(tt® tt)ii) = 2, 

((l/® Z/)l_l| (f, -T 2 +T^)\{V® l/)i_i> =2, 



((l®2)s|0?i-<7 2 f|(l®2)s) = (-85+12), for S = 0,1. 



(3.14) 



Hence, for spatially symmetric correlated pairs (i.e., T = 0, S = 1 and T = 1,S = 0) we have that in a simple 
schematic model with independent T = 1 and T = pairs 

S° T+ = (lu - wi) (^Af^l^o - 3 (^f(^)n + 7V^i-i)) . (3.15) 

This equation indicates that EWSR is influenced by both types of paring interactions, provided that the inequalities 
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L0 



^Wi and Af^T=o ^ 1 Jj\A" r )ii +AA (w) i-i) , (3.16) 



hold. 

If in Eq. (3.12) the places of operators Yj~ and Y~ are inter-changed, then one obtains the EWSR for the (3 + 
transition. It is easy to show that for the [3 + case, one has: 

S^~=S^ + . (3.17) 

Adding only an attractive T = force (luq > 0) to the mean field Hamiltonian, Ho, will increase 
Therefore, the EWSR for GT operators will increase while the NEWSR for GT operators will decrease. The same 
is true when we add only an attractive interaction in T = 1 channel. In this case it is interesting to note that the 
NEWSR decreases with increasing number of T — 1 nv pairs, while the EWSR does not depend directly on that 
number but increases with increasing number of proton-proton and neutron-neutron pairs. On the other hand when 
both T = and T = 1 attractive pairing forces are present, their effects tend to cancel in the EWSR for Gamow- Teller. 

IV. APPLICATION TO 44 TI WITHIN F 7/2 SHELL 

Here we want to evaluate numerically the contribution of various interactions on the GT strengths in a schematic 
single j— shell formalism. 

In the single j— shell model for the N = Z nucleus ^Ti (j = $7/2) the general form of a state with total angular 
momentum I, is: 

yal = £><"(./! J 2 ) [(j 2 ) Jl (f) J2 Y , (4-1) 

Ji,Ji 

where J\ is the total angular momentum of the two valence protons and J2 of the two valence neutrons. 

A. Non-energy weighted strength 

For the transition from a given 7 7r = + state to a given I v = 1+ state, the GT strength is: 



B+ T (a^(3) = 8\(j\\cr\\j)\' 
x V2J+ 1 



^D a0 {JJ)D^{JJ) 



J 



1 3 3 
3 J J 



(4.2) 



1 3 3 
3 J J 



If we sum over all final 1 + states, we get 

B+ T (a) = (8 - 6^,2) |0W>| 2 E \ Da °( JJ )\ 2 ( 2J + 1) 

j 

= 0.97959 [D Q0 (22)] 2 + 3.26531 [D a0 (44)] 2 + 6.85714 [L> Q °(66)f . (4.3) 

Here, for / 7/2 , | (j\ \a\ \j) | 2 = 72/7. 

Note that in the single j— shell model space there are four I n = + states, three with isospin T = and one with 
isospin T = 2. There are three I n = 1 + states, all with T = 1. We will now calculate Bq T for the J 71 " = + ground 
state of 44 Ti to 7 7r = 1 + T = 1 states in 44 V. We consider several cases whose B^ T values are given in Table 1. 

Case I. The isopairing Hamiltonian. 

We consider a simple isospin dependent Hamiltonian of the form 

H = -l + bh.t 2 , (4.4) 

which is equivalent to the Hamiltonian considered in Section 3 but written in the monopole form of French (see also 
Ref. [26]). When b > 0, the state r = 0+ T = 2 is pushed up in energy but still the three states F = 0+ T = are 
degenerate. In particular, when b = 4 we meet the situation of pure isopairing (T — pairing). 
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We again average over initial states 

B ir = \ E J2 B gtHT = 0)^/3). (4.5) 

a(T=0) /3 

In this case the total strength is obtained from the above equation by excluding from the summation over a the 
state a — \I^ = + ,T = 2). In order to do that we need to know the coefficient D T=2 ' J=0 (JJ), in the expansion 
(4.3). It turns out that these expansion coefficients are equal to the two particle fractional parentage: 

|D r =^=°(J,J)| = |((j 2 )J(i 2 )J|}/0)|. (4.6) 

(L> T=2 (0,0) = -0.5, D T=2 (2,2) = 0.3127, £> T=2 (4,4) = 0.5, L> T=2 (6,6) = 0.6009). This is due to the fact that the 
jtt = o+. T = 2 state in 44 Ti, can be obtained from the unique I* = + , T = 2 state in 44 Ca by applying on it twice 
the raising isospin operator. Hence, one finds, 

b £t = \ E s+^Cc - ^) = § I01k!b>l 2 C 1 - (0' 2 )-^0" 2 )^l>i 4 o) 

a(T=O),0 ' J 

x(2J+l){j J j J} =2.56. (4.7) 
Case II. Same as Case I but with b — 0. 

If we set b = then all four J = + states are degenerate. To the total strength in Case I we add the contributions 
of the J = 0+ T = 2 to all final J = 1+ T = 1 states. The added value is 

2|(i|H|j)| 2 X:(2J+l)|((i 2 )J(i 2 )J|}i 4 0)|(2J+l)|] ] j j}\ (4.8) 

When we average over the four initial states, we get 2.13. 

Now we consider another three cases obtained with pairing Hamiltonians which are different from the one given by 
Eq. (4.4). 

Case III. Standard (J = 0, T = 1) pairing. 

The case of isospin conserving ( J = 0, T = 1) pairing provides a non-degenerate ground state of 44 Ti by means of 
a schematic many body Hamiltonian fixed by the condition 

({f)'T\V\(f)i) = -5j, 5t,i MeV. (4.9) 

This is just the reduced isospin model of Flowers [27]. With condition (4.9), the matrix elements in the space of two 
protons plus two neutrons, which is the case of Ti, become: 





3 


J p] 




3 






3 


j p 


Hi 


3 




lo 





J 










H Jp ,j, = -25 Jp , 5j,, a - 4(2J P + l)(2j; + 1)1 j j J p W j j J'\. (4.10) 

loo o J lo J 

Taking into account the expressions of the 9-j symbols with one vanishing row one arrives at the following matrix 
representation of the many body Hamiltonian. 



H = -±- 
16 



33 VE 3 Vl3' 

V5 5 3V5 

3 3V5 9 3^ 

.a/13 V65 3Vl3 13 . 



(4.11) 



This Hamiltonian has a non-degenerate ground state, one excited state at 0.75 MeV and two degenerate states at 
2.25 MeV. The first eigenvalue is described by the state: 

$ 0+ = 0.8660[00]° + 0.2152[22]° + 0.2887[44]° + 0.3469[66]°. (4.12) 

The total strength for (3~ transition from the non-degenerate ground state + to any of the dipole states 1 + in the 
neighboring odd-odd nucleus 44 V, is given by Eq. (4.3) with D a0 (JJ) given in Eq. (4.12). The numerical value 
obtained can be seen in Table 1. 
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In the most often used case of like particle pairing and no irv pairing 

((j 2 n J \v\(fn J ) - {{fn J \v\(fn J ) = -sj, Mev, 

which violates isospin conservation, one has that the coefficients in Eq. (4.1) satisfy 

D°(J, J) = 6j t o 



(4.13) 
(4.14) 



and therefore one sees from Eq. (4.3) that Bq T = 0. Though at first sight this result is surprising, the reason for that 
is that we are introducing a pairing Hamiltonian in a single j— shell. Including only a single j— shell prevents pairing 
correlations to develop. So, we have to be careful in drawing too many conclusions from this particular result. 
Case IV. (T = 0, J = 1) pairing. 

In this case we consider a schematic many body Hamiltonian fixed by the condition 



((f)T\V\(f)l) = -5j, 1 6 Tfi MeV. 
For the J = + states of 44 Ti, the matrix elements become: 

Hj P ,j; = -f(J P )f{J' p ) , 

with 



f(J P ) 



i 2 ) Jp (f) Jp 



( j j Jp 

6(2J p + lWj j J p 



(4.15) 
(4.16) 

(4.17) 



I 1 1 

Numerically, the Hamiltonian matrix for J = + states in 44 Ti for this case is 



H = 



-0.18750 -0.33940 -0.20536 +0.22535' 

-0.33940 -0.61437 -0.37173 +0.40791 

-0.20536 -0.37173 -0.22491 +0.24681 

+0.22535 +0.40791 +0.24681 -0.27083. 



(4.18) 



Proceeding as in the previous case, we find after diagonalization Bq T — 2.46, which is very close to the number 
obtained in Case I where T — also dominate. 

The eigenfunctions and eigenvalues of this separable hamiltonian are easier to obtain than for Case III ( J = 0, T = 
1). The eigenvalue equation is 



-/ {Jp) E / (4) D K 4) = XD ( J p' J P ) ■ 

If we multiply by D(J p ,J p ) and sum, we get 



(4.19) 



Y,D{J P ,J P ) 



\J2D 2 (Jp,Jp) = A. 



It is easy to show that 

From the properties of 9j and 6j symbols, we find 

D ( J p , J p ) = Af 



f(Jp) 



(4.20) 



(4.21) 







/ 0.380 \ 


( 5/V5 ^ 




0.688 


| 69 




0.416 


\-2/yi3/ 




V -0.457/ 



(4.22) 



and then we get Bq T = 2.46. 

Case V. Equal (J = 0, T = I) and (J = 1, T = 0) pairing 
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((j 2 )t\V\U 2 )t) = - St.,o + 6j fi 5 T ,i) MeV. (4.23) 
For the J — states of 44 Ti the square matrix to be diagonalized is 



H 



-2.250 -0.47916 -0.39286 

-0.47916 -0.92687 -0.79100 -0.09598 

-0.39286 -0.79100 -0.78741 -0.42923 

-0.09598 -0.42923 -1.08333, 



(4.24) 



Note that the direct coupling ([0, 0] V[6, 6]) vanishes in this case. For (J = 0, T = 1) pairing this matrix elements is 
-0.22535 whilst for (J = 1, T = 0) pairing it is +0.22535. 
The ground state wave function is 

* = 0.8256[0, 0] + 0.4047[2, 2] + 0.3726[4, 4] + 0.1259[6, 6] . (4.25) 

We shall see that this is much closer to the realistic wave function to be discussed in the next section. The summed 
B(GT) strength is in this case 0.73. This is smaller than the values for pure (J = 0, T = 1) pairing and for pure 
(J = 1, T = 0) pairing. There is clearly a correlation of strength when both interactions are present and with the 
same sign. 

Case VI. MBZ 

This case corresponds to the situation considered by McCullcn et al. in [28] , where the ground state is also no longer 
degenerate. The quoted paper considers a more realistic two-body interaction whose matrix elements were fixed so 
that the spectrum in 42 Sc is reproduced. In MBZ one has the integrated effects of ( J = 0, T = 1) and (J = 1, T = 0) 
pairing, as well as a quadrupole-quadrupole interaction. For this interaction the non-degenerate ground state of 44 Ti 
is: 

* 0+ = 0.7608[00]° + 0.6090[22]° + 0.2093[44]° + 0.0812[66]°. (4.26) 

In the daughter odd-odd nucleus, the lowest level energy is a doublet I* = l + ,7 + at 0.6 MeV and a triplet V = 
2+, 3 + , 5 + at 1.5 MeV. The GT transition takes place from the ground state of 44 Ti, which is + , to the lowest state 
in 44 V, i.e. 1+, which has the structure: 

* 1+ = 0.9156I22] 1 + 0.3914[44] 4 + 0.0927[66] 1 . (4.27) 

Eqs. (4.26) and (4.27) determine the expansion coefficients D aI , which arc to be introduced in the expression (4.2) 
in order to calculate the Bq T strength for the transition + — > 1 + . 

Comparing the expressions (4.12) and (4.26), one sees that the coefficient Z?°(66) yielded in the schematic model 
described above, is much larger than that obtained by McCullen et al. [28]. The reason is that in ref. [28], the 
formalism involves implicitly the Q.Q interaction whose effect is to push up the state 6 + which results in diminishing 
the overlap coefficient I? (66). 

Summarizing, Case I corresponds to the situation of isopairing. Case II is associated to the situation of no in- 
teraction, i.e., all four J = + states are degenerate. Case III is the case of isospin conserving standard pairing, or 
Flowers model [27]. Case IV is the case of ( J = 1, T = 0) pairing, Case V is the case of an equal ( J = 0, T = 1) and 
(J = 1, T = 0) pairing and finally Case VI is the realistic calculation from ref. [28]. The results for the (3~ strengths 
corresponding to the cases described above, are collected in Table 1. 

We note that although the total strength for the 'no interaction' case is larger than that of T — isopairing, the 
averaged values are ordered differently. The reason is that the number of states for the first case (4) is larger than for 
the second case (3). The ordering of the total strengths given in the first column of Table 1, is consistent with the 
fact that the T = pairing decreases the total strength. In fact, increasing the number of T = pairs the blocking 
effect is increased and the occupation probability decreased, which results in a quenching effect for the total strength. 

A more drastic reduction is obtained with isospin conserving standard pairing (Case III) . An even stronger reduction 
is finally obtained in the realistic case considered in Case VI, which is caused by the fact that Q.Q interaction, implicitly 
involved in the case VI, favors approaching the limiting SU (4) symmetry for which Bq T would vanish. 

Note that for all cases considered, [3 + and f3~ strengths are equal. Hence, pairing interaction reduces both (3 + 
and j3~ strengths. These results complete those of ref. [17] showing that the two-body dipole interaction in the 
particle-particle channel suppresses part of the (3 + GT transition strength. Indeed, here we point out that a severe 
compression is taking place also for (3~ transition. 



11 



B. Energy weighted strengths 



We now consider results for several interactions: 
Case a. (J = 0, T = 1) pairing interaction 



{fy*) J T V {f^y) = -125(0)1*7,0 St,i ■ (4.28) 



The energy shifts and B{GT)' s for states in 44 Ti are as follows 

E(0+) = -2.25|£(0)|, 

E(l+) = -0.75|^(0)|, B(GT) Q+ ^ lt = 1.1507 
2?(l+) = 0, B(GT) o+ ^+=0 
£(!+) = 0, B(GT) Q+ ^ lt =0 



(4.29) 



EWS = 1.7260|£(0)|. 

Note that there is no B(GT) strength to the non-collective states, i.e., no strength to the states that are not shifted 
from their unperturbed position by the (J = + , T — 1) pairing interaction. 
Case b. (J = 1, T = 0) pairing interaction 

(fy2) J T V (fi /2 ) J T ) = -|25(l)|*/,i St,o. (4.30) 



Now we have 



E(0 + ) = -1.2976|£(1)|, 

E(l+) = — 1.02621^(1)1, BiGT)^^ = 0.03127 

E{1+) = -0.3190|£(1)|, B{GT) 0+ ^ 1 + = 2.4435 

£7(1+) = 0, B(GT) Q+ ^ lt = 



(4.31) 



EWS = 2.3996|£(1)|. 

Note that here also there is no strength to the non-collective state. In contrast to case a) ( J = 0, T = 1) pairing, 
we now have most of the strength going to the second 1+ state, not the first. This conclusion is consistent with 
our statement in Section 3 concerning the migration of the transition strength when T = and T = 1 pairings are 
considered separately. 

Case c. Equal (J = 0, T = 1) and (J = 1, T = 0) pairing 

f V 2) J T V (/ 7 2 / 2 )^) = -|2?(0, 1)| (<Wt,i + 5 Jtl S T . ) . (4.32) 



In this case 



£(0+) = -2.6622|£(0,1)|, 

E{1+) = -1.7131|£(0, 1)|, BiGT)^^ = 0.604 

E{1+) = -0.3826|£(0, 1)|, B{GT) Q+ ^+ = 0.122 

E(l+) = 0, B(GT) 0+ ^ lt = 



(4.33) 



EWS = 0.8523\E(0,1)\. 

This result confirms the statement made at the end of Section 3 saying that the T = and T = 1 pairings interfere 
destructively for EWS. 
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According to Eq. (3.15), S^ T = for ujq = ui\. However, here the value of EWS is small but not zero. The reason 
is that this case is different from the one considered in Section 3. 
Case d. MBZ 

E(l+) = 5.8146 McV, S(GT) 0+ ^ 1 + = 0.5180 
E(l+) = 8.5438 MeV, B{GT) Q+ ^+ = 0.0337 
£(1+) = 11.1393 McV, B(GT) 0+ ^ 1 + - 

(4.34) 

EWS = 3.2995 MeV. 

Note that results for the overall B(GT) strength for MBZ as well as the (J = + , T = 0) 44 Ti ground state wave 
function arc much closer to the results of the case of combined equal (J = + , T = 1) and (J — 1+, T — 0) pairing 
than they are to the individual pairings. Hence, one needs both types of pairing present in order to get realistic 
results. 

Note in particular that for (J = 1+, T = 0) pairing the [6,6] component of the (J = 0, T = 0) ground state 
wave function has the opposite phase of the same component in the more realistic MBZ wave function. The overlap 
< $(MBZ),$(J = 1+, T = pairing) > is only 0.76, whilst the overlap of MBZ with the mixed pairing case is 0.96. 



C. A note on the number of pairs 

For a system of n nucleons in a single j shell with total isospin T, we have the following results for the number of 
pairs, 

Total number of pairs: n(n — l)/2. 

Number of pairs with isospin T = 0: n 2 /8 + n/4 - T(T + l)/2. 
Number of pairs with isospin T = 1: 3n 2 /8 - 3n/4 + T(T + 1) /2. 

Hence, for the T = ground state of 44 Ti (n = 4) we have three T = pairs and three T = 1 pairs. 
Given a ground state wave function of the form 

* = ]T D(J P , J p ) [(il) J ° (jl) Jp ] J= ° (4.35) 

we can further obtain the number of pairs with total angular momentum J i2 (Ji 2 = 0, 1, 2, 3, 4, 5, 6, 7). The expression 
for the number of J i2 pairs is 



2|£(Ji2,Ji 2 )| 2 + |/(Ji 2 )| 2 (4.36) 



where 



/(J12) - 2(2 J 12 + 1) D(J P , J P )(2J P + 1)1 j j J P (4.37) 
Jp I J12 J12 J 



Note that D{ J 12 , J 12 ) is non-zero only for even (Ji 2 , T 12 = 1) states. 

We present results for the number of pairs with (J i2 = 0,Xi 2 = 1), (J i2 = l,Ti 2 = 0), and (J i2 = 7,Xi 2 = 0) in 
Table 2 for various interactions. The total number of pairs is of course 6, three with Ti 2 = and three with Xi 2 = 1. 

With the (J = 0, T = 1) pairing interaction we get 2.25 (Ji 2 = 0,Ti 2 — 1) pairs. Recall that the total number of 
T12 = 1 pairs is three. The fact that this interaction conserves isospin, forces the number of pairs to be less than 
three. There are only 0.25 ( Ji 2 = 1, T\ 2 = 0) pairs. The other 2.75 are shared by Ji 2 = 3, 5, 7. 

With a (J = 1,T = 0) pairing interaction the number of (J i2 = l,7i 2 = 0) pairs increases considerably to 1.297, 
whilst the number of (Ji 2 = 0,Ti 2 = 1) pairs decreases to only 0.433. When we have equal (J = 0, T = 1) and 
(J = 1, T = 0) pairing, we get 2.045 and 0.618 for the number of (J 12 = 0,Ti 2 = 1) and (J 12 = 1,T 12 = 0) pairs, 
respectively, result close to the MBZ result of 1.736 and 0.746. 

In general, there are more (Ji 2 = 7, T i2 = 0) pairs than (J i2 = l,Ti 2 = 0) pairs. The number of (J i2 = 7, T i2 = 0) 
is comparable to the number of (J 12 = 0,T 12 = 1) pairs. 

It is interesting to note that for the (J = 1, T = 0) pairing interaction, the wave function given by Eq. (4.22) can 
be written in another way, 
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l - P ni n 2 ) [(pmi) 



(p 2 n 2 ) Jl2=1 ] 



J=0 



(4.38) 



Jl2 = l 



where P is the exchange operator. In other words, first couple one neutron-proton pair to Ji 2 = 1 and then a second 
one. Couple the total to angular momentum zero and then antisymmetrize. At first, this seems like a very naive thing 
to do, but it works. Without the antisymmetry this looks like two J i2 = 1 pairs, but the act of antisymmetrization 
complicates things and gives us the result of Eq. (4.22). We get only 1.297 J12 pairs. It should be noted that for the 
(J = 0, T = 1) pairing interaction one cannot get a corresponding simple expression. 

V. THE EFFECT OF THE T=0 INTERACTION ON THE GAMOW-TELLER TRANSITIONS IN 44 TI 

Since the spin-flip configurations are important for Gamow- Teller transitions, inclusion of the f 5 / 2 orbit is essential 
to discuss on Gamow- Teller transitions as well as the Ikeda sum rule. In Sect. 4 the Ikeda sum rule is satisfied due 
to the relations B + (GT) = B~{GT) and Sp T+ = Sp T in spite of the limitations of the configuration space. Also 
without / 5 / 2 , a condensate of T — pairs cannot appear [25] . In this section we discuss the effect of T = interaction 
including the full fp shell. 

We give results for 44 Ti in which up to t nucleons are excited from the / 7 / 2 shell into higher shells. The t = case 
corresponds to a single j shell calculation and t = 4 to a full fp shell calculation. 
To find the effect of T = pairing we consider three interactions 

• II. The full FPD6 interaction [29]. 

• 12. Set the two particle matrix elements with isospin T = to zero but keep the T = 1 matrix elements of 
FPD6 unchanged. 

• 13. Set the T = matrix elements to -1.0 MeV but keep the T = 1 matrix elements of FPD6 unchanged. 

The idea is to see the effect of T = pairing by turning the T = interaction off and on. The 13 interaction has 
nonzero T = matrix elements but does not single out any particular ones. 

Just to be clear we distinguish between a constant matrix element and a constant interaction. A constant interaction 
a + bti ■ t 2 will have no effect on the wave functions or energy levels of a given total isospin. It will cause shifts in 
energies of states of differing isospin. A two-body matrix element can be written as 



If we have a constant interaction then this will vanish if (j 3 , j'4) is not the same as (ji, j 2 ). What we are doing however 
in 13 is setting all T = two-body matrix elements as above to -1.0 MeV. 

Those interactions have also been considered in Ref. [30]. The results are given in Table 3. From Table 3 one can 
see that in the single j shell calculation, the results for 12 and 13 are identical. This can be understood by the fact 
that in this small model space a constant matrix element and a constant interaction are identical. 

There is a large increase in B(GT) when we go from t = to t = 1. This is because we can now have transitions 
from f-j/2 to / 5 / 2 (there were not present in the t = calculation). 

But the greatest interest lies in comparing the results for II and 12. When the T = interaction is turned off (12) 
the values of B{GT) increase considerably, especially in the full calculation (t = 4). The change is from 1.236 to 
3.316. Putting it the other way, reintroducing the T = matrix elements (II) causes a large quenching of the summed 
B(GT) strength. This is also true in the single j shell calculation but it is not as pronounced. The pairing aspects of 
the T = interaction in FPD6 can be seen by examining the spectrum of 42 Sc where the excitation energies of the 
J = if and J = 7f T = states are at 0.6 MeV but the J = 3^ and J = 5f T = are at 1.5 MeV. We thus have 
the extreme spin states, lowest— J and highest— J, lower than the states with intermediate angular momentum. 

To see that this separation is important we introduce the interaction 13 which has constant T = matrix elements 
(E = —1.0 MeV). We see that the results for 13 are similar to those of 12 -large values of B(GT). Hence, the fine 
details of the full interaction (II), i.e., the pairing aspects for both low and high spin, are important. They are needed 
to explain the strong quenching of B{GT). 



In summary, we have investigated within a schematic model how T = and T = 1 pairing correlations show up in 
the energy weighted and non-energy weighted sums of [3 + and (i~ strengths. We find that both EWSR and NEWSR 




(5.1) 



VI. CONCLUSIONS 
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depend differently on the number of T = and T = 1 pairs in the ground state and a scheme is derived to assess 
the relative importance of T = 1 and T = pairing from experimental strength distributions obtained from (n,p) 
and (p, n) reactions. Formulas are given to get insight on pair numbers out of /3~ and /3 + strengths. In particular, 
we have found that attractive T — interactions alone tend to decrease the total GT strength, while increasing the 
EWS strength. Likewise ■nv T = 1 pairing alone. 

The numerical results presented for schematic j— shell calculations in the j3~ GT decay of 44 Ti show that indeed 
T = and T = 1 pairing correlations tend to decrease the GT strengths. The / 7 / 2 formalism is restricted and the 
sum strengths are smaller than the total strengths obtained in a major N— shell. Yet, this limited shell space allows 
to make simple calculations that illustrate the effect of the various interactions. The effect of T = pairing is best 
illustrated by the results shown on calculations in the full fp shell. 

In the future it will be interesting to apply this scheme to analyze both experimental data and realistic calculations. 
Calculations with more conventional forms of the pairing interaction [25,31] are at present being considered. 

To close this section we would like to say a few more words about previous results concerning the proton-neutron 
pairing and the single and double beta decay processes. As we mentioned, it was earlier noticed [17] that the transition 
amplitudes for (3 + and double beta decay are strongly influenced by the two body interaction in the particle-particle 
channel. Analyzing the second quantization expressions for this interaction, we notice that it is nothing else but 
the (T = 0, J = 1) pairing interaction when one deals with Gamow- Teller transitions, and the (T = 1, J = 0) 
pairing interaction when Fermi like transitions are considered [24] . We would like to stress that each of the processes 
mentioned above (Fermi and Gamow- Teller, (3~ and /? + ) is influenced by both types of pairing. 
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Table 1. Summed strengths for 44 Ti within / 7 / 2 shell with various schematic interactions (see text). 



Case 



~7F — 

GT 

Total 



rT 

GT 

Averaged over 
initial states 



I. Isopairing 

II. No interaction 

III. Isospin conserving (J = 0, T = 1) pairing 

IV. (J — 1, T — 0) pairing 

V. Equal (J = 0, T = 1) and (J = 1, T = 0) pairing 

VI. MBZ 



7.674 
8.531 



2.56 
2.13 
1.14 
2.46 
0.73 
0.55 



Table 2. Number of pairs 


with J 12 = 0, 1 and 7 in a single 


j shell calculation for 44 Ti. Results 


are for interaction III: 


Isospin conserving (J = 0, 


T = 1) pairing; IV: (J = 1, T = 


= 0) pairing; V: Equal ( J = 0, T = 


1) and (J = 1, T = 0) 


pairing; and VI: MBZ. 








Interaction 




Number of J12 pairs 






(Jl2 = 0, T12 = 1) 


(J12 = 1, Ti 2 = 0) 


(J12 = 7, T12 = 0) 


III 


2.250 


0.250 


1.250 


IV 


0.433 


1.297 


1.311 


V 


2.045 


0.618 


1.654 


VI 


1.736 


0.746 


1.948 



Table 3. The values of B{GT) in 44 Ti for the three interactions (11,12,13) and up to t nucleons excited from the / 7 / 2 
shell. 



II 



12 



13 



0.627 
3.192 
1.771 
1.236 



0.840 
4.109 
3.442 
3.516 



0.840 
3.495 
3.122 
3.495 



17 



